Bubbles in ^^Si and 220? 



00 
O 

o 

(N 
D 

in 



> 

o 
00 
o 



M. Grasso,^ L. Gaudefroy,^ E. Khan,^ T. Niksic,'^ J. Piekarewicz,^ O. Sorlin,^ N. Van Giai/ and D. Vretenar^ 

' Institut de Physique Nucleaire, Universite Paris-Sud, IN2P3-CNRS, F-91406 Orsay Cedex, France 
^ CEA/DAM lie de France, Service de Physique Nucleaire, 
Bruyere-le-Chdtel, F-91297 Arpajon Cedex, Prance 
"'Physics Department, Faculty of Science, Zagreb University, 10000 Zagreb, Croatia 
'^Florida State University, Tallahassee, Florida 32306, USA 
^Grand Accelerateur National d'lons Lourds (GANIL), GEA/DSM-CNRS-IN2P3, 
Bd Henri Becquerel, BP 55027, F-l\o76 Caen Cedex 5, France 

Bubble nuclei are characterized by a depletion of their central density. Their existence is exam- 
ined within three different theoretical frameworks: the shell model as well as non-relativistic and 
relativistic microscopic mean-field approaches. We propose ^*Si and ^'^O as possible candidates for 
proton and neutron bubble nuclei, respectively. In the case of ^^O, we observe a significant model 
dependence, thereby calling into question the bubble structure of ^^O. In contrast, an overall agree- 
ment among the models is obtained for '^''Si. Indeed, all models predict a central proton density 
depletion of about 40%. This result provides strong evidence in favor of a proton bubble in '^''Si. 
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I. INTRODUCTION 



The "bubble" structure of atomic nuclei is character- 
ized by a depleted central density. Although it is some- 
what unexpected that a "hole" can be made in a nuclear 
system where nuclear forces generate a saturation density 
(Pq ~ 0.16 fm~'^), this phenomenon has been discussed 
for many decades. Indeed, the possibility of bubble nu- 
clei started with the pioneering work of Wilson in the 
40s [l| who studied the low-energy excitations of a thin 
spherical shell, up to the first microscopic calculations of 
Campi and Sprung in the 70s [1| . More recently, bubbles 
have been discussed in superheavy and hyperheavy nu- 
clei [3, 0] ■ The promise of producing more exotic nuclei 
with the new generation of RIB facilities has revived the 
interest in this subject. 

Due to the absence of a centrifugal barrier, s-orbitals 
have radial distributions peaked in the interior of the 
nucleus, with their corresponding wave function extend- 
ing further into the surface depending on the number of 
nodes. In contrast, orbitals with non-zero angular mo- 
menta are suppressed in the nuclear interior and do not 
contribute to the central density. Therefore, any vacancy 
of s-orbitals is expected to produce a depletion of the 
central density. By using electron scattering from ^'^^Pb 
and ^os^p^ large momentum transfers, the radial 

distribution of the 3s proton orbital was experimentally 
mapped and shown to closely resemble the one predicted 
by an independent particle model. The agreement ex- 
tends from the center of the ^°^Pb nucleus all the way to 
the surface and reproduces accurately the nodal struc- 
ture of the wave function Q . Differences in the charge 
density between ^°^Ph and ^°^T1 revealed that about 80% 
of the proton removal strength came from the 3s state, 
thereby leading to a significant depletion of the proton 
density in the nuclear interior. Specifically, the depletion 



fraction, defined as 



(1) 



amounts to F = 11(2)%. Note that here pc and pmax 
represent the values of the central and maximum charge 
density in ^°^T1, respectively. Yet the small energy dif- 
ference between the 3s 1/2 and the 2^3/2 proton orbitals 
plus the coupling of the 3s 1/2 proton to collective excita- 
tions in ^°^Pb, yield a proton hole strength in ^ostI that 
is shared among the 3s 1/2 and 2^3/2 orbitals, with the 
former carrying about 70% of the strength and the latter 
the remaining 30%. Consequently, the central depletion 
in ^oS'p^ relative to ^°^Pb is not as large as if the full hole 
strength would have been carried by the 3s orbital. Using 
similar arguments, the depletion at the center of ^°'*Hg 
is not expected to be very large, as the two-proton hole 
strength will be again shared among the 3s 1/2 and 2^3/2 
orbitals. Therefore, the search for the best bubble candi- 
dates should be oriented towards nuclei with an s orbital 
well separated from its nearby single-particle states and 
where correlations are weak. This latter feature arises 
mainly for nuclei located at major shell closures. 

Recently, the formation of a proton bubble due to 
the depletion of the 2si/2 orbital was investigated in 
^^Ar 0, Q and in the very neutron- rich Ar isotopes [1]. 
In *^Ar the proton 2si/2 and 1^3/2 orbitals are almost 
degenerate. Thus, as in the case of ^°^Pb, pairing corre- 
lations will lead to a significant occupancy of the 2si/2 
orbital [91], thus weakening the bubble effect. This weak- 
ening will continue to hold for any = 28 isotone be- 
tween Z = 20 and Z = lA as long as the 2si/2 and 
orbitals remain degenerate, as shown for instance 
in Fig. 3 of Ref. [l^l- For very neutron-rich Ar isotopes, 
such as ^^Ar, the si/2 proton orbital is predicted to move 
significantly above the state, hindering the role of 
pairing correlations @, |ll| . Unfortunately, the produc- 
tion of this exotic nucleus is far beyond the present and 
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near-future capabilities of RIB facilities. 

A more suitable region of the chart of the nuclides to 
search for a proton bubble is that of the = 20 iso- 
tones. Between Z = 20 and Z = 16 the Si/2 orbital is 
located about 6.5 MeV above the ^5/2 orbital and about 
2.5 MeV below the ^3/2 orbital, thereby forming two sub- 
shell closures at Z — 14 and Z — 16, respectively [12]. In 
addition, the iV = 20 shell closure is rigid enough to hin- 
der significant coupling to collective states. Assuming 
a sequential filling of proton orbitals, the 2si/2 orbital 
should be completely empty in '^^Si while fully filled in 
^^S. This may lead to an important change in the proton 
density distribution between '^^S and '^^Si making ■^^Si an 
excellent candidate for a bubble nucleus. Concomitantly, 
both Skyrme and Gogny Hartree-Fock-Bogoliubov mod- 
els predict a spherical shape for -^"^Si [H, [l3|- Other pos- 
sible candidates in the Si-isotopic chain, such as ^^Si and 
^^Si, are not optimal as they are deformed [H, [lB|- For 
these nuclei several correlations hinder the development 
of a bubble. The mirror system of the ("^^S, '^^Si) system, 
(■^^Ca, '^^Ca), could not be studied at present because the 
■^^Ca nucleus has so far not been observed. 

A neutron bubble may be found in the oxygen chain, 
where large A = 14 (between c?5/2 and Si/2) and 
A = 16 (between Si/2 and d^/'?) subshell gaps of about 
4.2 MeVflTl. [Til and 4 MeV [l^, respectively, have been 
determined. Combined to the large pro ton gap at Z = 8, 
the [13, US lU and ^^Q [il,[ii,l23 theferore behave 
as doubly magic nuclei. In this case the change in the oc- 
cupancy of the 2si/2 neutron orbital will occur between 
2^0 and 2^0, making 2^0 a good candidate for a neutron 
bubble nucleus. 

The present article aims at determining whether '^'^Si 
and 2^0 could be considered as good proton and neu- 
tron bubble nuclei, respectively. Various theoretical ap- 
proaches will be employed to test the robustness of the re- 
sults. In Sec. II these nuclei are analyzed in terms of shell- 
model calculations so that occupancies of the proton and 
neutron orbitals may be determined. In Sec. Ill we start 
by addressing (in III-A) the role of pairing correlations in 
mean field approaches, followed then with results on nu- 
cleon density profiles obtained from: (i) non-relativistic 
Hartree-Fock (HF) and Hartree-Fock-Bogoliubov (HFB) 
(in III-B) and (ii) relativistic mean field (RMF) and rela- 
tivistic Hartree-Bogoliubov (RHB) (in III-C) microscopic 
calculations. Comparison to experimental data will be 
made whenever possible. Conclusions are drawn in Sec. 
IV. 



II. SHELL-MODEL PREDICTIONS 

The occurrence of bubbles in nuclei, as previously de- 
fined, is directly linked to the occupancy of Si/2 orbitals. 
For both bubble candidates under study in this article, 
220 and ^^Si, experimental values for the occupancies 
are not yet available. Thus, we rely hereafter on Shell 
Model (SM) calculations that are known to give reli- 



TABLE I: Ground state occupation numbers of neutron or- 
bitals obtained in the SM framework for ^'*0 and ^^O. 



Orbital 


24q 


22q 




5.75 


5.38 




1.89 


0.34 




0.36 


0.28 



able estimates in neutron-rich nuclei occupying the sdpf 
shells 

Shell-model calculations have been performed with 
the ANTOINE code [H, ^ using the SDPF-NR inter- 
action [2^ . The full sd valence space was considered for 
both proton and neutron excitations for each nucleus un- 
der consideration. Moreover, in the calculations for '^^S 
and "^^Si, ApAh neutron excitations were allowed from the 
vd'^/2 to the vf-!/2 orbitals. 



A. Neutron Bubble: ^^O and ^^O 

The mean occupation numbers of neutron orbitals de- 
duced from SM calculations are reported in TableHl Con- 
trary to naive expectations, the mean number of neu- 
trons removed from the j^Sx/2 orbital while moving from 
2'*0 to 2^0 amounts to only 1.55 (or 78% in fraction of 
shell coccupancy). The remaining neutron strength is 
removed from the vd^/2 orbital (0.37 or 18%) and the 
i'd^/2 orbital (0.08 or 4%). These numbers suggest that 
the depletion of the neutron density in the nuclear in- 
terior of 2^0 relative to 2^0 may not be as large as 
required for the formation of a bubble structure. To 
quantify the size of the neutron hole the radial depen- 
dence of the neutron wave functions in 22,24q jjg^yg been 
calculated using a standard Woods-Saxon potential ^2^, 
the parameters of which are as follows: Vb — -50 MeV, 
a — 0.65 fm and tq = 1.25 fm. The neutron densities 
displayed in Fig. [1] are the sum of the squared wave func- 
tions calculated with the Woods-Saxon potential (with- 
out spin-orbit term) weighted by the occupation numbers 
obtained within the SM framework. The orbitals occu- 
pied by the core neutrons (namely, the \vsx/2, lt'P3/2 
and l^'j3i/2 orbitals) are assumed to be completely filled. 
The figure shows the variation of the neutron density be- 
tween 2^0 and 2^0. In particular, the predicted depletion 
fraction, as defined in Eq. ([1]), is found to be F = 25%. 
As will be shown in the following section, the depletion 
in 2^0 is considerably smaller than in •^"'Si. One notices 
that in the lighter Oxygen isotopes this depletion factor is 
not expected to increase as the neutron li'dc^/2 orbital is 
depleted in concert, thus lowering the maximum value of 
the density around the surface of the nucleus. It follows 
that the relative difference of the density in the vicinity 
of the surface and at the interior of the nucleus is also 
reduced. 
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TABLE II; Same as Table|T]for ^''S and ^*Si. The third column 
shows the experimental occupancies obtained in Ref. [27|] for 





Orbital 




^^Si 


^''^S (Ref. [27j) 


7rld5/2 


5.85 


5.75 


6.1(12) 


7r2si/2 


1.88 


0.09 


1.63(32) 


7rld3/2 


0.27 


0.16 


0.31(6) 



B. Proton Bubble: ^'^S and ^''Si 

The mean occupation numbers of the proton 1^3/2, 
2si/2 and 1^5/2 orbitals in ^^S have been obtained from 
the ^^S(d,'^iJe)^^P experiment done by S. Khan and col- 
laborators 122]. The sum of the deduced spectroscopic 
factors from the proton pickup reaction from the whole 
sd states amounts to Y^C^S « 7.9. Within the 20% 
uncertainties of the method, this exhausts the complete 
spectroscopic shell model strength of ^C^S* = 8. Ex- 
perimental occupancies of the 2si/2 and 1^3/2 orbitals 
as reported in Table HIl are 1.63 and 0.31, respectively. 
Note that the small occupancy of the 1^3/2 state is due 
to correlations. 

As already mentioned, ApAh neutron excitations from 
the to the vf'j/2 orbitals have been allowed in the 

present SM calculations. This provides an estimate of the 
contribution of neutron cross-shell excitations (across the 
iV = 20 shell closure). Such contributions are found to 
be negligible (< 2%). Indeed, both and ^''Si ground 
states have quasi pure single-particle wave functions. For 
34Si rs], the {ird,/2f{vsdY^ [{^d./^fijrsy^fii^sdY^] 
configuration represents about 87% [85%] of the ground 
state wave function. The mean occupation numbers for 
the proton orbitals deduced from the present SM calcu- 
lations are reported in Table |TT1 The agreement with the 
experimental values for ^®S is very good, lending confi- 




FIG. 1: Neutron density of ^"O and ^^O obtained usins the 
occupation numbers of neutron orbitals resulting from SM cal- 
culations and radial wave functions obtained with a Woods- 
Saxon potential (see text). 



dence to the SM predictions for ^ Si. The study of the 
'^^Si(d,'^i/e)'^'^Al reaction channel should be used in the 
future to determine if the occupancy of the 2si/2 proton 
orbital has indeed dropped to nearly zero in ^"^Si. 

The proton fraction removed from the 27rsi/2 orbital, 
while moving from ■^^S to ■^'^Si, is significantly larger than 
the corresponding case in the Oxygen isotopes. More- 
over, a larger mean occupation number of the c?5/2 or- 
bitals is predicted in ^"'Si as compared to ^^O. Both of 
these effects combined to explain the larger depletion in 
the nuclear interior of '^'*Si, as displayed in Fig. [51 The 
densities presented in this figure have been obtained in 
the same way as for ^*'^^0. The depletion factor in '^''Si 
is found to be F = 43%, significantly larger than in ^^O. 
Thus, we conclude that within this approach, '^'^Si is a 
good candidate for a proton bubble. We are aware that 
the nucleon densities reported above may not be as re- 
liable as the mean occupation numbers extracted from 
the SM code. The self-consistent microscopic treatment 
presented in the following sections is more appropriate to 
provide radial densities. 

III. MEAN FIELD CALCULATIONS 

Self-consistent mean-field approaches enable to deter- 
mine microscopically the density distributions of nuclei. 
Nucleon occupation factors may also be determined by 
taking pairing correlations into account. As a first step 
in describing the density distributions, one must deter- 
mine the role of pairing correlations, if any, on the devel- 
opment of proton and neutron bubbles in 34 Si and 22 Q 
nuclei, respectively. 

A. Pairing effects 

As already alluded to in the Introduction, ^^O is ex- 
pected to behave almost as a doubly-magic nucleus, be- 
ing that the = 14 sub-shell closure has been experi- 
mentally determined. However, as shown in the previous 

1 0x1 0-2^ 




FIG. 2: Same as Fig.[T]for proton densities in ^'^S and ^"'Si. 
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section, SM calculations predict a 17% occupancy of the 
2s neutron state, suggesting that pairing correlations are 
likely to have some effect on this nucleus. The effect of 
pairing correlations on the neutron density profile of ^^O 
will be shown in the following subsections for both the 
non-relativistic and the relativistic mean-field cases. 

Let us now consider the case of '^''Si. To start, we 
discuss the role of pairing in the non-relativistic case. 
In this case pairing correlations can be modeled in the 
Skyrme-Hartree-Fock-Bogoliubov (Skyrme-HFB) frame- 
work by adopting the following zero-range density depen- 
dent pairing interaction: 



1 



p(r) 
Po 



S (ri - r2) ; 



(2) 



with r] = 0.5 (mixed surface- volume interaction), a ~ 1 
and Po = 0.16 fm~'^. In the particle- hole channel, we 
employ the SLy4 Skyrme parametrization that is well 
suited to describe neutron-rich nuclei. 

We fix the parameter Vb in Eq. (2) to reproduce the 
two-proton separation energy in '^'^Si. Note that the two- 
proton separation energy is defined as: 



S2p = E{N, Z)-E{N,Z^2), 



(3) 



where E{N, Z) is the total binding energy of the {N, Z) 
nucleus. It should be noted that the experimental value 
of 33.74 MeV is already reasonably well reproduced with- 
out pairing: the HF value is equal to 35.19 MeV. More- 
over the HFB calculations — which include the pairing 
interaction — yield negligible corrections, a,s Z = 14 is 
predicted by the HFB approach to be a robust sub-shell 
closure in agreement with the shell model spectroscopic 
factors (see Table II where the SM occupation of the s 
state is only 4.5%). Thus, we can safely perform the anal- 
ysis of this nucleus with the HF model, where pairing is 
absent. Note that the above results will be confirmed in 
Subsec. C within the relativistic approach where pairing 
effects are also found to collapse in the case of ^''Si. 
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FIG. 3: HF neutron densities (in units of fm"^) of ^^O (full 
line) and (dashed line) calculated with the Skyrme in- 
teraction SLy4. The dotted line represents the SLy4-HFB 
neutron density of ^^O. 



in Fig. 3). Note that the density profile of ^^O remains 
unchanged when pairing is switched on. 




FIG. 4: HF proton densities (in units of fm"^) 01^*^8 (dashed 
line) and "^''Si (solid line) calculated with the Skyrme interac- 
tion SLy4. 



B. Non-relativistic mean field approach 

In this Subsection both •^■^O and "^"'Si are analyzed as 
possible candidates for bubbles nuclei. Figure 3 displays 
neutron density profiles in ^^O (full line) and ^"^O (dashed 
line) calculated self-consistently within the SLy4-HF ap- 
proach. In the case of ^''O the neutron single-particle 
state 2si/2 is assumed to be fully occupied. The deple- 
tion of the central density in ^^O relative to ^"'O is clearly 
visible. However, the bubble profile is not evident: since 
the central neutron density in ^^O is strongly enhanced, 
the depletion in ^^O does not lead to the development of 
a significant central hole. The central depletion fraction 

is ~ 13%, much weaker than the SM result. As one 
switches on pairing and chooses the same parameters as 
in Ref. [l^ for the pairing interaction, the central hole 
is seen to be completely washed out (see the dotted line 



The SLy4-HF proton density profiles calculated in •^"'Si 
and '^^S (where the s state is fully occupied) are shown 
in Fig. 4. One observes that the bubble is much more 
prominent in this case than in ^^O. The depletion frac- 
tion is ~ 38% and in very good agreement with the 
SM value. We should mention that pairing is expected to 
modify the density profile of '^^S, but this is not relevant 
here since this density is plotted only to better appreci- 
ate the central hole in ■^^Si. By comparing the HF proton 
density in '^*Si with the HF neutron density in ■^^O, one 
observes that the central value in •^"'Si is much lower than 
in ^^O. The contribution to the central value of the den- 
sity is entirely due to the first s wave function, i.e. the 
Is. The difference between the two central values may 
be related to the presence of a neutron excess at the sur- 
face of '^'*Si. The effect of this neutron-skin on the proton 
Is 1/2 wave function is to attract and push it towards the 
surface, thereby lowering its value at the center. This 



5 



can be observed in Fig. 5 where the neutron (proton) Is 
contribution to the HF density is plotted for ^^O ("^"'Si). 



0,08 



0,06- 




FIG. 5: Neutron (proton) Is contributions to the density (in 
units of fm-3) for ^^O (^''Si). 



C. Relativistic mean field approach 

As in the previous section, calculations are performed 
for the two Oxygen isotopes ^^O and ^^O as well as for 
the two A'' = 20 isotones ^''Si and '^^S, but now using a 
relativistic mean field (RMF) approach. Pairing effects 
are evaluated within the relativistic Hartree-Bogoliubov 
(RHB) model. In one particular realization of the rel- 
ativistic formalism the dynamics of the system is dic- 
tated by an interacting Lagrangian density of the follow- 
ing form: 



9s4>- {. 



A 



5p, 



^+2(1 + ^3)^^ 



2 



r 



(4) 



where ?/> represents an isodoublet nucleon field interact- 
ing via the exchange of two isoscalar mesons — a scalar 
(0) and a vector (F^) — one isovector meson (&^), and 
the photon (A'^) [23, [SOl- In addition to meson- nucleon 
interactions, the Lagrangian density is supplemented by 
non-linear meson interactions with coupling constants de- 
noted by K, A, and Av that are responsible for a soft- 
ening of the nuclear-matter equation of state, both for 
symmetric and pure-neutron matter. For the RMF case 
we consider two parametrizations: the very successful 
NL3 parameter set [3l|, H^l and a more recent set known 
as FSUGold [S^ . The main difference between these two 
models lies in the prediction of the density dependence 
of the symmetry energy. This difference manifests it- 
self in significantly larger neutron skins for NL3 than for 
FSUGold [3^. Neutron skins for the two isotones of in- 
terest in the present work, alongside other ground-state 
properties, have been listed in Table Hill for ^''Si and ^^S. 



Model 


B/A(MeV) 


i?ch(fm) 


7?„-i?p(fm) 


NL3 


8.36 


3.13 


0.25 


FSUGold 


8.28 


3.13 


0.21 


Experiment 


8.34 






NL3 


8.50 


3.26 


0.12 


FSUGold 


8.42 


3.26 


0.09 


Experiment 


8.58 


3.28 





TABLE III; Binding energy per nucleon, charge radii, and 
neutron skin thickness for ^*Si (upper block) and "^^S (lower 
block) as predicted by the two RMF models used in this work. 
When available, experimental data is provided for compari- 
son. 
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FIG. 6: RMF neutron densities of ^^O and ^''O calculated 
with the two RMF models described in the text. 



RMF neutron densities for the two neutron-rich iso- 
topes ^^O and ^''O are displayed in Fig. 6. Whereas 
the RMF results show a mild model dependence, differ- 
ences between the relativistic and non-relativistic models 
are significant. Indeed, in contrast to the non-relativistic 
case, the relativistic results display no enhancement of 
the central neutron density in ■^'^O. Moreover, the removal 
of both 2si/2 neutrons from ^^O yields a strong deple- 
tion of the interior neutron density in ^^O. As a result, a 
central depletion fraction of F = 34%(28%) is predicted 
for by the FSUGold(NL3) model. These values are 
significantly larger than the 13% depletion fraction ob- 
tained with the SLy4-HF parametrization, but close to 
the SM expectations. It may be interesting to elucidate 
the source behind this discrepancy which might be re- 
lated to the saturation mechanism in RMF and HF cal- 
culations. 

In the case of the larger nuclei "^''Si and "^^S one ob- 
serves, now in agreement with the non-relativistic and the 
SM cases, how the proton density of "^''Si is significantly 
depleted in the nuclear interior and how the proton bub- 
ble disappears as soon as the 2si/2 proton orbital is filled 
in ^^S (see Fig. [7]). Indeed, not only is the central density 
in 34si strongly depleted, but the central density in 
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FIG. 7: RMF proton densities of ^"^S and ^*Si calculated with 
the two RMF models described in the text. 



FIG. 8: Neutron density profiles of ^^O and calculated in 
the RHB model with the density-dependent RMF interaction 
DD-ME2, and Gogny pairing. 



also appears to be get greatly enhanced. This behavior 
results in central depletion factors for ^''Si of F = 42% 
and F — 37% for the FSUGold and NL3 parameter sets, 
respectively. These numbers are in good agreement with 
the non-relativistic prediction of F 38% and SM of ^ 
43%. 

Let us quantify now the effects of pairing correlations 
within the RHB model. A medium dependence for a rel- 
ativistic mean-field interaction can either be introduced 
by including non-linear meson self-interaction terms in 
the Lagrangian, as in the case of NL3 and FSUGold, 
or by assuming an explicit density dependence for the 
meson- nucleon couplings. This is the case of the DD- 
ME2 model [33| that we adopt to perform RHB calcula- 
tions. The couplings of the a-meson and w-meson to the 
nucleon are assumed to be of the form: 



9iip) = 9i{Psat)h{x) for i = CT,w 



where 



l + b,{x + dif 



(5) 



(6) 



is a function of x = p/ Psat, and psat denotes the nucleon 
density at saturation in symmetric nuclear matter. Con- 
straints at nuclear matter saturation density and at zero 
density are used to reduce the number of independent 
parameters in Eq. ([6]) to three. Three additional param- 
eters in the isoscalar channel are gaiPsat), gui{Psat), and 
•nia — the mass of the phenomenological a meson. For the 
p meson coupling the functional form of the density de- 
pendence is suggested by Dirac-Brueckner calculations of 
asymmetric nuclear matter: 



9p{p) = 9p{Psat) exp[-ap(x - 1)] 



(7) 



and the isovector channel is parametrized by gpipsat) and 
a p. Bare values are used for the masses of the lo and p 



mesons: rriu 



783 MeV and = 763 MeV. DD-ME2 



is determined by eight independent parameters, adjusted 
to the properties of symmetric and asymmetric nuclear 
matter, binding energies, charge radii, and neutron radii 
of spherical nuclei [3J]. The interaction has been tested 
in the calculation of ground state properties of large set 
of spherical and deformed nuclei. When used in the rel- 
ativistic RPA, DD-ME2 reproduces with high accuracy 
data on isoscalar and isovector collective excitations [3J? 
]• 

In Figs. [8] and [9] we display the neutron (proton) den- 
sity profiles for 22,24o^ 34gi 36g^ calculated in RHB 
model [35i with the DD-ME2 effective interaction in the 
particle-hole channel, and with the Gogny interaction j36j 
in the pairing channel 

1=1,2 

- H,P^ - AhP^P^), (8) 

with the set DIS [sj for the parameters /i^, Wi, Hi, 
and {i = 1,2). 

For and '^^Si the RHB calculation with the DD- 
ME2 interaction predicts neutron and proton density pro- 
files similar to those calculated with NL3 and FSUG- 
old. Because of the large gaps: between vsi/2 and 1^^3/2 
(> 4 MeV) in ^"^O, and between 1^(1^/2 and 7rsi/2 (> 6 
MeV) in '^^Si, we find a pairing collapse in these nuclei, 
in agreement with non-relativistic predictions. On the 
other hand, the inclusion of pairing correlations has a 
pronounced effect on the neutron and proton density pro- 
files in ^^O and ^^S, respectively. When pairing is set to 
zero (dash-dot curves in Figs. [8] and [9|) the or- 
bital is empty in ^^O, and the 7rsi/2 orbital is fully occu- 
pied in ^^S. The resulting DD-ME2 density profiles are 
again very similar to those calculated with the two other 
RMF interactions. However, the pairing interaction in 
the RHB model calculation modifies the occupancy of the 
two 2si/2 orbitals, thus reducing the pronounced bubble 
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FIG. 9: Proton densities of and ^''Si calculated in the 
RHB model with the DD-ME2 interaction plus Gogny DIS 
pairing. 



Nucleus 


SM 


SLy4 


SLy4 


NL3 


FSUGold 


DDME2 


DDME2 






HF 


HFB 


RMF 


RMF 


RMF 


RHB 


22q 


25% 


13% 





28% 


34% 


29% 


10% 


34Si 


43% 


38% 


38% 


37% 


42% 


36% 


36% 



TABLE IV: Central fraction of depletion F 



in the neutron density of O, as well as the prominent 
cusp in the proton density of "^^S. 

In the DD-ME2 model the F values are found equal to 
29%, 10% and 36% for ^^O (without pairing), 22 Q (with 
pairing) and ^*Si (same result with and without pairing), 
respectively. 
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